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Abstract 


Work during the period July 1 - December 31, 1985, has concentrated 
on the application of the equations derived in the preceding period to 
the maneuvering and vibration suppression of the Spacecraft Control 
Laboratory Experiment (SCOLE) model. Two different situations have been 
considered: 1) a space environment and 2) a laboratory environment. 

This report covers the first case and consists of a paper entitled 
"Maneuvering and Vibration Control of Flexible Spacecraft", presented at 
the Workshop on Structural Dynamics and Control Interaction of Flexible 
Structures , Marshall Space Flight Center, Huntsville, AL, April 22-24, 
1986. The second case will be covered in the report for the next 
period. 
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ABSTRACT 


This paper is concerned with slewing a large structure in space and 
suppressing any vibration at the same time. The structure is assumed to 
undergo large rigid-body motions and small elastic deformations. A 
perturbation method permits a maneuver strategy independent of the 
vibration control. Optimal control and pole placement techniques, 
formulated to include actuator dynamics, are used to suppress the 
vibration during maneuver. The theory is illustrated by simultaneous 
maneuvering and vibration control of the Spacecraft Control Laboratory 
Experiment (SCOLE) model in a space environment. 

I. INTRODUCTION 

The problem of simultaneous maneuver and vibration suppression of 
spacecraft is becoming increasingly important. Some projected NASA 
missions involve experiments consisting of the control of flexible 
bodies carried by a shuttle in an earth orbit. Other missions involve 
laboratory simulations of similar experiments. The equations of motion 
for both types of experiments have been presented previously (Ref. I). 

A perturbation technique permitting the maneuver strategy to be 
formulated independently of the vibration problem was also presented. 

In a following investigation, a straightforward rotational maneuver 
strategy and an efficient vibration simulation technique were presented 
(Ref. 2). Rotational maneuvers were shown to cause structural 
vibrations. The purpose of this paper is to develop ways of suppressing 
these vibrations during the maneuver by means of feedback control. 

Turner and Junkins (Ref. 3), Turner and Chun (Ref. 4) and Breakwell 
(Ref. 5) addressed the problem of rotational maneuvering and 
simultaneous vibration suppression of flexible spacecraft for two- 



dimensional models. In all cases, the methods used represent extensions 
of rigid-body maneuvering techniques, requiring the solution of a two 
point boundary-value problem. The maneuver and vibration control 
problems are coupled and numerical difficulties arise as the order of 
the system increases (Ref. 5). Baruh and Silverberg first suggested 
separating the maneuver and vibration control problems (Ref. 6). 

However, the vibration control did not include feedback on the rigid- 
body modes, so that the spacecraft orientation could not be corrected 
during the open- loop maneuver. 

The equations of motion of a flexible spacecraft consist of six 
nonlinear ordinary differential equations for the rigid-body motion of a 
reference frame attached to the spacecraft in undeformed state coupled 
with a set of partial differential equations for the vibration of the 
elastic members relative to the rigid frame. Hence, the equations 
describing the motion of a flexible spacecraft during a certain maneuver 
represent a set of nonlinear hybrid differential equations. In general, 
hybrid systems of equations do not permit closed-form solution, so that 
one must consider an approximate solution, which implies spatial 
discretization and truncation. 

The nonlinear equations of motion can be solved by a perturbation 
approach. The approach consists of separating the equations into a set 
of equations for the rigid-body motions, representing zero-order 
effects, and a set of equations for the small elastic motions and 
deviations from the rigid-body motions, representing first-order 
effects. The perturbation technique permits a maneuver strategy that is 
independent of the elastic vibration. 


The order of the perturbation equations for the vibration 
suppression is often so large that some reduction is necessary. To this 
end, the elastic motion can be expanded into a series consisting of 
premaneuver eigenvectors acting as admissible vectors. These admissible 
vectors clearly do not decouple the equations of motion, but using a 
smaller number than the order of the system permits a reduction in the 
order. We refer to the corresponding equations as quasi-modal 
equations . The task of simulation and vibration control can be carried 
out conveniently by means of the quasi-modal equations. 

At each sampling time, the control forces are formulated as if the 
premaneuver eigenfunctions are the true eigenfuctions at that time. 

This results in a modeling error, but this error does not create a 
stability problem for robust control techniques, such as natural 
control, which do not require the exact eigenfunctions in control 
formulation (Refs. 7-11). 

Any maneuver can be regarded as a single-axis rotation, where in 
general the axis of rotation is not a principal axis. The single-axis 
maneuver has the advantage of simplifying the matrices of time-dependent 
coefficients in the perturbation equations. The rigid-body maneuver 
represents an open-loop control strategy for minimum-time, single-axis 
rotation about an axis which is not necessarily a principal axis. 

The vibration control is carried out in discrete time, which 
amounts to regarding the system as having constant coefficient over the 
duration of any sampling period. This permits the use of control 
algorithms developed for time-invariant systems, such as optimal control 
and pole placement. Actuator dynamics can degrade the performace of the 
control system. The inclusion of actuator dynamics requires a 



reformulation of the vibration control techniques mentioned above. The 
theory is illustrated by simultaneous maneuvering and vibration control 
of the SCOLE model in a space environment. 

II. EQUATIONS OF MOTION 

We consider a space system consisting of a shuttle carrying an 
antenna connected to the shuttle by means of a mast, as shown in Fig. 

1. The shuttle is assumed to be rigid and the mast and antenna are 
deformable. The motion of the spacecraft is referred to a given 
reference frame XQygZQ embedded in the rigid shuttle (Fig. 1). The 
reference frame has six degrees of freedom, three rigid-body rotations 
and three rigid-body translations. We propose to derive the equations 
of motion by means of the Lagrangian approach, which requires the 
expressions for the kinetic energy, potential energy and virtual work. 

Considering Fig. 1, the position of a point S in the rigid shuttle 

relative to the inertial frame XYZ is R^ = R + r. Moreover, the position 
of a point A on the elastic appendage is R^ = R + a + u, where u is the 
elastic displacement. The velocities of these points are then 

Rg = R + u) x r (1) 

R. = R + u) x (a + u) + u (2) 

~M ~ ~ 

where R and u are the translational and angular velocities of the frame 
x 0 y 0 z 0 w ^ th res P ec t to the inertial frame, respectively, and u is the 
elastic velocity of the point relative to XQy 0 z 0 . Hence, the kinetic 
energy of the spacecraft is 

T-jJ |R| 2 du s +ij lR A | Z dm A (3) 

"S *A 

where m<j and m A are the masses of the shuttle and appendage. 


respectively. 



The potential energy is due to the combined effects of gravity and 


strain energy. The gravitational potential can be expresssed as 
V = - Gm [J | R + r | _1 dm Q + J |R + a + u| _1 dm fl ] 


where m e is the mass of the_earth and G is the universal gravitational 
constant. The strain energy can be expressed as an energy inner product 
denoted by [ , ] (Ref. 12). This includes centrifugal and gravitational 
stiffening effects on the appendage. The total potential energy then 


becomes 

V = |[u,u] + V g (5) 
Denoting by f<- and f A the force vectors per unit volume of the shuttle 
and appendage, respectively, we can express the virtual work as 


6W = J f $ -6R s dD $ + J f A -6R A dD A (6) 

D S °A 

where and D A are the domains of the shuttle and appendage, 
respectively. The system is discretized in space by expressing the 
elastic displacements vector in the form 

u = $q (7) 

where $ is a matrix of space-dependent admissible functions and q is a 
vector of time-dependent generalized coordinates. The kinetic and 
potential energies and virtual work can be expressed in matrix form by 
introducing a rotational transformation matrix C from the XYZ frame to 
the *g.yo 2 0 f r ame, where the elements of C are nonlinear functions of a 
set of Euler angles a^, a^, a^. Lagrange's equations of motion can then 
be written in the symbolic form 


_d 

f— 1 

+ iV - 

c t f 

dt 

aft 

sR 
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f — 1 

u + 

8 V 

dt 

l~J 

3a 

8a 
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(8a) 


(8b) 



where the matrix D is defined by the expression 


ui = 0 (a)a (9) 

The resulting equations of motion are nonlinear due to the large rigid- 
body motion. 


III. PERTURBATION METHOD 


Consider a first-order perturbation on the quantities R and a 
R = Rg + R^ , a = ciq + (10a, b) 

where the first-order terms R^ and are small compared to the zero- 
order terms Rq and ctg. Introducing Eqs. (10) into the nonlinear 
equations of motion and separating orders of magnitude, we obtain zero- 
order and first-order perturbation equations. Before proceeding with 
this technique, we first develop some expressions relating the pertur- 
bations a^, a^, ^ in the Euler angles, a^, a^, with sma ^ angular 

deflections 6^, B^, 8^ expressed in the body-fixed frame. This will 
permit expressing all the variables in the perturbation equations in 
terms of components along the Xq^O axes. 

It can be shown that the vector of body-fixed perturbation angles 
can be related to the vector of perturbed Euler angles by 


8 = D(ctg)a^ (11) 

The perturbed angular velocity vector can then be expressed as 


~T. 


u> = u)q + - ojgS + B 


(12,13) 


In keeping with our objective of expressing the first-order perturbation 
equations in the body-fixed frame, Eq. (10a) is replaced by 


3 = 5o + c o?i 


(14) 


where R^ is now a vector measured with respect to axes XgygZg and the 



rotation matrix has the perturbed form 

C = Cq + = bCq (15a, b) 

The applied forces and moments can also be expressed in first-order 
perturbed form as follows: 

F = F q + Fi , M = M Q + M x (16a, b) 

The zero-order equations of motion, which govern the structure as 
if it were rigid, can be expressed as 


+ C 0 S 0~0 + C 0 W 0 S 0~0 + 


Gm 


73 l m Bo + (' - 3 ?o?J> c J§ol 


~T ' ^ m e ~T • ~T 

S 0 C 0~0 + ,3 S 0 C 0?0 + *0-0 + “O^-O = -0 


(17a) 

(17b) 


where Rq is a unit vector in the direction of Rq, Iq is the mass moment 
of inertia matrix about point 0, m is the mass of the spacecraft and 


§0 


- v 


dm c + f 
S J m, 


a dm. 


(18) 


The first-order linear perturbation equations, which govern the 
vibrational motion of the structure, can be expressed in the matrix form 
Mx + Gx + (K $ + K ns )x = F* (19) 

where 


; T - i*I 


T T, f *T IF T 
6 q ] , F = l F 1 


9o + 9i> 


(20a, b) 


M = 




m 0 0 
0 m 0 
0 0 
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(20c, d) 
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M. = r $ T « dm. , j(w) = L (aw + [au])« dm. 

A Tn^ A ~ in ^ ~ M 


IV. RIGID-BODY MANEUVER 


(20e) 


(20f) 


(20g) 


(21a, b) 
(21c, d) 
(21e,f) 


The perturbation method permits the maneuver strategy to be 
designed independently of the vibration control. Indeed, a straight- 
forward single-axis minimum-time maneuver strategy is developed in Ref. 
2. It is shown in Ref. 2 that the axis of rotation need not be a 
principal axis, so that any general rotational maneuver is possible. 

The maneuver policy is formulated in continuous time but implemented in 
discrete time, so that some error can occur. Also, the equations 


governing the maneuver are nonlinear, so that the solution of Ref. 2 can 
be unstable. However, the vibration control includes feedback control 
of the rigid-body modes and can stabilize the spacecraft as well as 
reduce the error caused by discrete-time sampling. 

The most desirable control technique for a maneuver excites only 
the desired rigid-body motion, not the elastic modes. From Ref. 2, the 
components of the maneuver force distribution exhibiting these 
characteristics can be expressed as 

F^p.t) = x(p)m(p)e 2 (t) (22a) 

F 2 (p,t) = -z(p)m(p)e(t) ' (22b) 

F 3 (P.t) = y(p)m(p)e(t) (22c) 

where e(t) is the desired angular motion, m(p) is the mass density at 
point p, and x(p), y(p) and z(p) are the components of the position 
vector of point p with respect to the center of rotation. The actuating 
forces are proportional to rotational rigid-body modes, so that they 
will not excite the elastic modes and cause undesirable vibration. Of 
course, distributed forces can only be implemented approximately with 
discrete actuators, which tends to excite some vibration. Also, 
centrifugal forces can cause vibration during the maneuver, so that 
vibration control may be necessary. 

V. QUASI-MODAL EQUATIONS OF MOTION 
During the maneuver, the gyroscopic and stiffness matrices, and 
hence the eigenvalue problem, are functions of time. In Ref. 2, a 
truncated set of the premaneuver eigenvectors is used as a set of 
admissible vectors to simplify and reduce the order of the equations of 
motion to a form called the quasi-modal equations of motion. Using this 
approach, Eq. (19) can be reduced to the quasi-modal form 


(23) 


u(t) + G(t)u(t) + [a + K(t)]u(t) = f(t) 

where 

G(t) = X T G(t)X, K(t) = X T K t (t)X (24a, b) 

are reduced-order gyroscopic and stiffness matrices, u is a vector of 
the quasi-modal coordinates defined by 

x ( t) = Xu(t) (24c) 

where X is a rectangular matrix of the lower premaneuver eigenvectors 
normalized so that 

X T MX = I , X T K Q X = A (24d ,e) 

in which Kq contains the constant terms in the stiffness matrix and a is 
a diagonal matrix of the premaneuver eigenvalues, and 

f(t) = X T F*(t) (24f) 

is the vector of modal forces. Comparing Eq. (23) with Eq. (19), we see 
that the premaneuver eigenvectors have not decoupled the equations of 
motion. However, as the maneuver velocity decreases, the time-varying 
terms decrease in magnitude and the equations approach an uncoupled 
form. Also note that the mass matrix has been reduced to the identity 
matrix, which is convenient for casting Eq. (23) in state space form. 

VI. VIBRATION CONTROL 

The tangential and centrifugal disturbance forces during the 
maneuver can cause vibration to be suppressed by means of feedback 
control, as demonstrated in Ref. 2. The perturbation method permits the 
vibration control to be formulated separately from the maneuver control, 
resulting in the quasi-modal equations of motion, 
a) Actuator Dynamics 

The equation of motion for the rth mode can be expressed as 
u r (t) + ^u r (t) = f dr (t) * f r (t) 


( 25 ) 



where f^ r (t) is the sum of the rth modal disturbance and maneuver 
control forces and f f (t) is the rth modal vibrational control force 
defined as the rth component of the vector 

f(t) = X T F ( t ) (26) 

where F(t) is a vector of control forces and moments. The modal 
disturbance force f dr (t) includes terms arising from G(t) and K(t), as 
can be concluded from Eq. (23). These terms are negligible when the 
maneuver angular velocity is small relative to the lowest natural 
frequency of vibration and the maneuver angular acceleration is 
similarly small. The equation governing the output of a typical set of 
actuators is assumed to have the vector form 

F(t) = aF(t) + bF (t) (27) 

~ - ~ L 

where a and b are constants and F ft) is the command force vector. 

~c ' 

Multiplying Eq. (27) by X T , the equation for the rth modal control force 
can be expressed as 


f r (t) . »f r (t) + bf cr (t) (28) 

where f r and f cr are the rth components of the vectors 

f = X T F (29a) 

and x 

f c = X T [ C (29b) 


respectively. In discrete time, if the feedback force f r is computed at 
t.j, the earliest it can be applied is the next sampling time t i+1 . 

However, considering Eq. (28) in discrete time, f cr is the quantity that 
is updated at t i+ j and the desired force f r cannot actually be applied 
until t^. This problem can be circumvented by including the actuator 
dynamics, Eq. (27), in the control formulation. 

Using Eq. (25) and assuming that f dr (t) can be neglected, we can write 

( 30 ) 


o . 

U + u) U = f 

r r r r 



Then, introducing Eqs. (25) and (28) into Eg. (30), the combined dynamic 


equation can be expressed as 

?. o 

u = -a) u + au + aw u + bf 
r r r r r r cr 

Equation (31) can be rewritten in the state form 


z = A z + bf 
~r r~r ~ cr 


where 


i} - [u u u 1 
~r 1 r r r J 


A = 
r 


0 

0 

r 

at/ 


1 

0 


b T = [0 0 b] 
0 
1 
a 


(31) 

(32) 
(33a, b) 

(33c) 


r r 
b) Optimal Control 

The minimum-time optimal control policy is bang-bang and involves 
determining three-dimensional switching surfaces for each mode. At each 
sampling time the modal state is to be located with respect to the 
switching surfaces, and based on this the command modal force is set at 
the positive or negative maximum value. The calculation of three- 
dimensional switching surfaces is not a simple matter and using them in 
the above mentioned fashion can be computationally time consuming (Ref. 
13). If the measured modal states are noisy, the modal forces could be 
switched frequently. This implies frequent, abrupt structural 
accelerations which, when applied with spatially-discrete forces, tend 
to destabilize the uncontrolled modes. 

As an alternative, we consider a quadratic performance measure in 
conjunction with the independent modal-space control (IMSC) method. In 
particular, we consider a performance index consisting of a weighted sum 
of the elastic energy and the control effort (Refs. 8, 9). This 
performance index has the form 



( 34 ) 


J = z J 
r=l r 


where 


J = J (z^Q z + 
r j q v ~r N r~r 


R f )dt 
r cr' 


(35) 


are modal performance indices, in which 


Q 


r 



0 0 
1 0 
0 0 


(36) 


is a weighting matrix and R r is a weighting factor. For the elastic 

? 

modes, q r is taken as the rth eigenvalue and for the rigid-body 
modes q r is chosen on the basis of pole-placement considerations. As 
the weight R r is decreased, the rate of modal energy dissipation is 
increased and, of course, more effort is required. Hence, R f can be 
chosen on the basis of the available control command force for vibration 
suppression. The command forces required by this performance index tend 
to be smooth functions. Smooth forcing time histories are preferable in 
that they do not tend to excite the uncontrolled modes to the extent 
that forces containing discontinuities do. 

To minimize J, each modal performance index J r can be minimized 
independently. The steady-state Riccati equation for the rth mode is an 
algebraic matrix equation of order three and can be expressed as 


A - 

A T K - 

- (D + ■ 

J- K bb T K 

• r 

r r 

M r 

R r r~~ r 

i symmetric 

matrix 

to be determined; 

k rl 


symm. 


k r2 

V4 



k r3 

k c 
r5 

k c 
r6 



(37) 


( 38 ) 



The linear, state feedback control law is of the form 


f = - i- b T K z (39) 

cr R r - r~r v 

Introducing Eqs. (33a, b) into Eq. (39), the control law becomes 

f = - Jr- (k _u + k c u + k c u ) ' (40) 

cr R v r3 r r5 r r6 r' K 

r 

which minimizes the performance index given by Eq. (35) (Ref. 13). In 
solving Eq. (37), the required entries of K r can be shown to satisfy the 
following expressions: 


k , = 
r3 

-S (auy + d ) 

b z r r 

(41a) 

k 4 ,C. 
r 6 4 

+ k^ c C-> + k 2 c C 0 + k ,C, + C n = 0 
r 6 3 r 6 2 r 6 1 0 

(41b) 

k c = 
r5 

h 2 9 

-ak r + -§ 5 — k , 
r 6 2 R^ r 6 

r 

(41c) 

where 

d r 

la 2 u) r + q r b 2 /R r l ^ » c q = J" ( ak r 3 “ 2 ^ 

r 

(42a, b) 

C 1 = 

1 fA j. a 2\ r b 2 , 2 . 2. 

- (d + a«o ), C 2 - 2d R (a> r + a ) 

(42c, d) 

C 3 

r r r 

.4 .6 

-b a r _ b 

? 9 3 

OA D C H Q A D° 

(42e,f ) 


are constant coefficients which can be evaluated for each mode off- 
line. The control law given by Eq. (40) is IMSC modified to include 
actuator dynamics. IMSC is also called natural control because the 
closed-loop modes are identical to the open-loop or natural modes, so 
that natural coordinates remain natural after control (Ref. 9). Natural 
control is efficient because it wastes no effort in changing the shape 



of the modes (Ref. 14), as other control techniques do. In this regard, 
we recall that only eigenvalues affect stability and not eigenvectors. 
Equation (40) can be expressed in the form 

f = -g ,u - g 0 u - g o u (43) 

cr 3 rl r 3 r2 r 3 r3 r v ' 

where the control gains are defined as follows: 


bk 


r3 


bk 


r5 


J rl 


r2 


J r3 


bk 

~R 


r 6 


(44a-c) 


r r r 

Introducing Eq. (43) into Eq. (32), considering Eqs. (33) and expanding 
the characteristic determinant, we conclude that the closed-loop poles 


must satisfy the equation 

s 3 + (bg r3 - a)s 2 + (bg f2 + u, 2 )s + (bg rl - aJ) = 0 


(45) 


c) Pole Allocation 

In optimal control, a performance measure is defined, perhaps 
arbitrarily, and minimized. The pole-allocation method is a modal 
control method in which the poles are chosen for each mode, again 
somewhat arbitrarily, and the actuator forces are computed to produce 
these preselected poles. The most general form of the poles for the rth 
mode defined by Eqs. (32) and (33), and whose characteristic equation is 
represented by Eq. (45), can be expressed as 


s rl ' a r + 1 B r’ S r2 a r 1B r’ S r3 Y r 


(46a-c) 


where and y^ are related to the time constants for the rth mode 
and e is the closed-loop modal frequency. The characteristic equation 
associated with the poles given by Eqs. (46) is 

S 3 + (-y r - 2a r )s 2 + (2y r a r + a 2 + S 2 )s + (-y^ 2 - y^ 2 ) = 0 (47) 
Comparing Eqs. (47) and (45), we obtain the control gains 

g rl = 1 law r " Y r (a r + B l )] 

g r2 = 1 (2Y r a r + 4 + B r ' * g r3 = E (a “ 2a r ' Y r } 


(48a) 
(48b, c) 



Hence, the poles given by Eqs. (46) can be chosen for each mode and 
implemented with the command modal force given by Eq. (43) and the gains 
of Eqs. (48). 

d) Implementation of Modal Control 
i . Modal actuation 

Implementation of any IMSC technique requires an inversion of Eq. 
(29b), so that the applied command force F c can be found from knowledge 
of the modal command force f . For discrete actuators, the simplest 
implementation technique is a projection method (Ref. 9). Premulti- 
plying Eq. (24d) by the pseudo-inverse of xJ t denoted by X - ^, we obtain 
X- T = MX (49) 

Premultiplying Eq. (29b) by X'"'” and considering Eq. (49), we obtain 

F c = MXf c (50) 

which is the discrete counterpart of a distributed generalized force and 
must be applied with actuators at every finite element node. This may 
require an impractically large number of actuators, so that only an 
approximate "distributed" force can be applied. In the projection 
technique, distributed control is implemented approximately by means of 
discrete actuators (Ref. 9). 

Another method of calculating the applied command force F c requires 
a numerical inversion (Ref. 8). The advantage of this technique is that 
it produces a subset of the desired components of the modal force vector 
exactly, rather than only approximately as in the case of projected 
control. Hence, this type of control is a true natural control in that 
the closed-loop eigenfunctions of the controlled modes are identical to 
the open-loop eigenfunctions. Note that the inversion must be 
recalculated in the event of an actuator failure. The forces computed 



from both the inverse and projected methods tend to excite the 
uncontrolled modes giving rise to so-called control spillover (Ref. 8, 

9). 

i i . Modal estimation 

Modal state estimation is another common problem inherent in the 
implementation of modal control techniques. This problem is somewhat 
analogous to that of modal force implementation. It requires the 
estimation of u(t) from measurement of x(t), which requires premulti- 
plication of Eq. (24c) by the matrix product X^M, so that 

u(t) = X T Mx(t) (51) 

which is the second half of the expansion theorem (Ref. 12). Equation 
(51) and its time derivatives form what are known as modal filters (Ref. 
15). Application of modal filters presupposes either "distributed" 
sensors or sensors at every node for a finite element model. An 
acceptable approach to modal state estimation is to interpolate the 
discrete sensor measurements by means of admissible functions to produce 
approximate displacement and velocity profiles, where the latter are 
functions of the spatial variables. These functions are then introduced 
into the modal filter equations to yield the approximate modal states. 

An alternative approach to modal state estimation involves a 
numerical inversion of a matrix. In a manner analogous to the inversion 
technique for force implementation, the displacement vector and modal 
matrix of Eq. (24c) can be partitioned as follows: 




where x is an s x 1 vector of measured displacements, x is the 
~s ~ns 

complement of x $ , B $ is an s x c submatrix of the modal matrix X and B 
is the complement of B_. Hence, the equation 


ns 


x = B u 
~s s~ 


(53) 


relates the modal displacements with the measured displacements. 

Premultiplying Eq. (52) by B“^, we obtain 

u = B' 1 ^ (54) 

where B~^ is only a pseudo-inverse of B s if s * c. The expression 

x = B B -1 x (55) 

~ns ns s ~s v ’ 

can be regarded as observation spillover. Equation (55) implies 
assignment of values to the unmeasured components of the displacement 
vector. 


e) Output Feedback Control 

Modal control techniques have become quite common in the field of 
control of structures because of their ability to take advantage of the 
physical and mathematical properties of the natural modes of vibration. 
The main drawback common to all modal control methods is the problem 
encountered in modal state estimation with discrete sensors. Direct 
feedback, where the command control force is related directly to the 
measured state through control gains, avoids modal estimation 
entirely. In addition, if this feedback force is the physical (rather 
than modal) force, then the modal force implementation is also 
avoided. Hence, a useful control technique might take advantage of the 
concepts of natural modes but be applicable to output feedback. To this 
end, the relationship between linear, output feedback control and linear 
modal feedback control will now be explored. 



Consider a linear control law having the special form 

F c (t) = -A[x(t)g 1 + x(t)g 2 + x(t)g 3 ] (56) 

where g^, g? and g 3 are control gain factors, F c is a vector of actual 
command forces, x, x and x vectors of measured displacements, velocities 
and accelerations and A is a weighting matrix of constant coefficients. 
Using Eq. (24c) we obtain 

F c (t) = -AX[u(t)g 1 + u(t)g 2 + u(t)g 3 ] (67) 

From Eq. (43), however, we conclude that, if the modal control gains are 
the same for all the modes, then the modal control force can be 
expressed in the vector form 

f Q (t) = -u(t)g 1 - u(t)g 2 - u(t)g 3 (58) 

Hence, Eq. (57) takes the form 

F c = AXf c (59) 

Recalling that Eq. (50) also relates f and F , and Eq. (50) is true 
regardless of the control technique used as it is an integral part of 
the expansion theorem, it is obvious that the weighting matrix A should 
be set equal to the mass matrix. Hence, if A is replaced by M in Eq. 
(56), the output feedback control law takes the form 

F c (t) = -M[x(t)g 1 + x(t)g 2 + xg 3 l (60) 

so that the control specified by Eq. (60) is equivalent to the modal 
control of Eq. (58). However, the control gain factors remain to be 
determined. 

A method of determining a particular set of control gains having 
good physical and mathematical basis is known as uniform damping control 
(Ref. 16). To this end, we consider Eqs. (46), which represent the 
desired closed-loop poles for pole allocation. But, changing the 
natural frequencies can be regarded as a waste of effort in the case of 



vibration suppression. In view of this, we choose s r in Eqs. (46) as 
being equal to the natural frequency w of the rth mode. The expansion 
theorem states that the motion of the structure can be represented as a 
linear combination of all the modes. Hence, it is reasonable to force 
all the modes to decay at the same rate of time, which implies that 
cy = a and v r = y in Eqs. (46). Introducing these values into Eqs. 

(48), the gains for the rth mode can be expressed as 
9 r I = l (a - x)uy - yd 2 ], g r2 = £ (2yq + a 2 ), g r3 = (a - 2a - y) 

(61a-c) 

If y is set equal to the decay rate of the actuator response a, then the 
modal gains for all the modes become 

g^ = -aa 2 /b, g 2 = (2 aa + a 2 )/b, g 3 = -2a/b (62a-c) 

which are identical for all the modes. Introducing Eqs. (62) into Eq. 
(60), we obtain a linear feedback control law providing uniform decay 
rate for all the modes without altering the natural frequencies. Hence, 
modal estimation and implementation is bypassed entirely, although the 
formulation takes advantage of the concepts of modal control. In fact, 
uniform damping control can also be derived as a first-order 
approximation of natural control if 

R . _b£_ 

2a 2 « 2 (63) 

and a « to for all modes, 
r 

As with other distributed control techniques, uniform damping 
control can only be implemented approximately with discrete actuators. 
However, this control method has the advantage of being applicable in a 
decentralized and collocated sense, and one that is known to be robust 
(Ref. 16). If an actuator and sensor pair is located at node i, acting 



in direction 2 , the control command force at that location takes the 
form 


F . = -m. (x. a, + x. q„ + x. q_) 

Cl 2. 1 2 ' 12^1 U a 2 U y 3 ; 


( 64 ) 


■where m. is an entry of the mass matrix, and x. , x. and x. are 

1 2, l x* 1 it i x 

correspondents of the displacement, velocity and acceleration vectors, 
respectively, all of which correspond to node i and direction 2 . This 
is known as decentralized control because the force at a point is 
related only to measurements taken at that point and, of course, is only 
possible when the actuators and sensors are collocated. Comparing Eqs. 
(60) and (64), we conclude that the approximation involved in 
decentralized control consists of ignoring the off-diagonal terms in the 
mass matrix. If collocation is not possible, the measured state must be 
interpolated spatially to approximate the state at the actuator 
locations. Another advantage of uniform damping control is that all the 
modes are controlled, even with discrete actuation, not just a subset as 
with typical modal control techniques. 

The mass matrix has been found to be the weighting matrix both for 
the open-loop maneuver forces, Eqs. (22), and for a particular closed- 
loop vibration control force, Eq. (60). There is a certain physical 
content in this result. Obviously, the uniform acceleration of a 
spacecraft demands more force at more massive sections. For the 
maneuver, this acceleration rate is the desired rotational acceleration 
and for the vibration control, it is the required vibrational decay 
rate a. 


Mil. NUMERICAL RESULTS 


The SCOLE configuration of Fig. 2 was modeled by means of the 
finite element method. The nodal locations, each representing six 



degrees of freedom, are denoted by the small circles. The mast 
supporting the antenna is a steel tube 10 feet long. The antenna 
consists of 12 aluminum tubes, each 2 feet long, welded together to form 
a hexagonal -shaped grid. The shuttle is simulated by a steel plate of 
uniform thickness with a mass of 13.85 slugs. The details of the model 
can be found in Ref. 2. The natural frequencies of the model in a space 
environment and cable detached are given in Table 1. 

The maneuver strategy of Ref. 2 was applied to the rigid-body model 
of the spacecraft. The differential equation governing the actuator 
behavior was assumed to be of first-order, as in Eq. (27), with a = - 10 
and b = 1. The histories of three rotational maneuvers about the Xq 
axis are presented. The discrete-time switching history of a 0° to 30° 
minimum-time rotation with M max = 20 Ib-ft is illustrated in Fig. 3a. 
Note that the acceleration overshoots the target state because of the 
discrete-time switching. Figures 3b and 3c illustrate the continuous- 
time switching histories of 0° to 180° rotations with M max = 20 lb-ft 
and M max = 60 lb-ft, respectively. 

To illustrate the effects of the actuator dynamics, the model is 
struck with an impulsive force of magnitude 5 lb in the yg direction. 

The structure is initially in a state of rest. The ensuing motion is 
modeled with 15 modes and the first 10 modes are controlled with natural 
control applied with distributed actuators and sensors. Both the 
classical form of natural control and natural control adapted to include 
the actuator dynamics is used to suppress the vibrations. 

Figures 4a and 4b contain plots of the total energy of the system 
versus the control performance for various values of the actuator 
response decay rate a. The values of R in the performance functional 



are 0.5 and 0.1 for Figs. 4a and 4b, respectively, while q = 1 in both 
cases, for all modes. The symbol + signifies a solution with classical 
natural control, in which actuator dynamics is not included. In all 
cases, the performance is improved with inclusion of actuator dynamics 
in the control formulation. Of course, as a -*• - ® this difference 
approaches zero. Note that, comparing Figs. 4a and 4b, this difference 
also depends on the value of R. As R decreases, the modal decay rate 
increases according to Eq. (63). Hence, the benefits derived from 
including the actuator dynamics in the control formulation depend on the 
ratio of the desired modal decay rate to the actuators' decay rate. Of 
course, as this ratio is decreased, the effect of actuator dynamics 
becomes less important. As can be observed from a comparison of Figs. 
3a-c, a similar conclusion can be made about the effect of actuator 
dynamics on the maneuver control. In this case, the ratio of interest 
is that of the mean angular maneuver velocity to the actuator decay 
rate. 

The remainder of this section is concerned with the vibration 
control of the 30° and 180° roll maneuvers in space. The first fifteen 
modes are modeled for the simulations. Simulations are presented 
without vibration controls in Ref. 2, showing the excitation of the 
structure by centrifugal and tangential forces. 

Figure 5 is a time-lapse plot of the spacecraft during the 30° roll 
maneuver without vibration control. The rotation is produced by 
actuators located on the shuttle, so that centrifugal and tangential 
forces cause structural vibration. The view is from directly behind the 
structure, showing the y 0 z 0 plane with the x Q axis directed into the 

At each plotting sampling time, two plots appear, one in dashed 


paper. 



lines representing the structure as if it were rigid and the other 
representing the deformed structure. As the structure is accelerated, 
the appendage lags behind its desired position and then it bounces 
forward to precede the desired configuration during deceleration. When 
the maneuver ends, the appendage continues to vibrate about the desired 
target state. 

The control techniques of Sec. 6 were applied to control the 
vibration of the spacecraft during the 30° roll. In all the presented 
cases, R = 0.01 and q = 1.0 for all the controlled modes. The time 
constant of the actuators' equations of motion was assumed to be a = 

-10, as in the maneuver strategy. Figure 6 is a time-lapse plot of the 
spacecraft during the maneuver with uniform damping control suppressing 
the vibrations of the first 9 modes using the ten actuators, six on the 
shuttle and four thrusters, two at the end of the mast and two at the 
antenna hub. Comparing Figs. 5 and 6, we conclude that uniform damping 
control during the maneuver provides excellent performance. 

The spacecraft was maneuvered through the 30° angle with both 
distributed actuators and the 10 actuators discussed above. These 10 
actuators were used for both maneuver control and vibration suppression. 
First, natural control approximated with projected actuating and 
distributed sensing was used. Next, 10 sensors were collocated with the 
actuators and the maneuver was repeated using natural control 
approximated by projected sensing and actuating. Finally, uniform 
damping control was used with the same 10 collocated actuators and 
sensors. 

Figures 7a and 7b show the total energy in the system and the total 
command effort versus time, respectively, for all four of the above 


mentioned control strategies. At the completion of the maneuver, energy 
is lower in all four strategies and the performance during the maneuver 
is always better compared to the case of no vibration suppression shown 
in Ref. 2. Distributed sensors and actuators yield the best performance 
by far. Note that uniform damping control performs almost as well as 
projected actuating and distributed sensing. This is to be expected, as 
the former is an approximation of the latter. However, it must be 
recalled that uniform damping control is decentralized and requires only 
a finite number of sensors, equal to the number of actuators. 

The case of projected actuating and sensing is clearly inferior to 
uniform damping control, even though both techniques were applied with 
the same collocated sensors and actuators. The difference lies in the 
feedback method. Uniform damping control is decentralized, so that each 
sensor supplies information only for the collocated actuator, requiring 
no interpolation. On the other hand, the projected actuating and 
sensing strategy is centralized so that the modal states are estimated 
through interpolation of signals from all the sensors. The actuators 
are commanded from the estimated modal states. Hence, to reduce modal 
estimation error with projected sensing, more sensors are needed, as 
discussed in Sec. 6. Better modal filtering techniques can help to 
reduce the modal estimation error, but in general, more sensors than 
actuators are required. 

The main control objective of the SCOLE project is aiming the 
antenna within a certain tolerance in minimum time. Hence, a study of 
the antenna rotations versus time is of interest. Figures 8a-c 
illustrate the instantaneous antenna rotations with respect to the frame 
XqYqZq about the x Q , y Q and z Q axes, respectively. The plots include 



the responses produced during the maneuver with all four of the above 
control implementation techniques. The relative quality of performance 
is as mentioned above. 

The effect of controlling only a finite set of modes is reflected 
in Fig. 8b in the form of residual motion of the uncontrolled modes. 
Uniform damping control dissipates energy in all of the modes, not just 
a finite set. From Fig. 8b, however, we observe that it is not very 
effective in controlling certain modes due to actuator placement. 
Nevertheless, with time, it does remove all of the energy from the 
system, unlike other implementation techniques. 

In reality, because of saturation, actuators can only produce a 
finite amount of force or moment. The effect of limited command forces 
was investigated for the 30° roll maneuver and uniform damping 
control. The shuttle command forces and moments were limited to 2 lb 
and 20 ft-lb, respectively. The four thrusters on the appendage were 
each limited to 0.15 lb. These values correspond to saturation of the 
thruster on the antenna hub in the yg direction. Figure 10 contains 
plots of this command force versus time for limited and unlimited 
actuation. The corresponding plots of modal vibrational energy versus 
time in Fig. 11a show the adverse effect of command force truncation. 

The price for force truncation is increased vibrational energy as well 
as an associated increase in overall effort, as illustrated by Fig. lib, 
although the control scheme is still effective. 

The maneuver and control techniques are now demonstrated for a more 
extreme situation, a 180° roll maneuver with M max = 60 ft-lb. From Fig. 
3c, we conclude that this is a relatively severe maneuver in that a 180° 
roll, starting and ending at rest, takes place in less than 4 seconds. 



Attempting this maneuver with actuators located only on the shuttle and 
without vibration controls resulted in deformations of the appendage 
exceeding greatly the small elastic motions assumption. 

Natural control and uniform damping control were employed to sup- 
press the vibration of the structure during the maneuver using R = .01 
and g = 1.0. First, natural control was applied with distributed 
actuators and sensors to establish a "best case". Uniform damping 
control was then applied with the same 10 actuators and sensors used 
previously. Finally, uniform damping control was used with two 
additional thrusters (and collocated sensors), one in the Zg direction 
at the hub of the antenna and the other at the extreme tip of the 
antenna in the yg direction. 

The time-lapse plots of the structure during the maneuver showed no 
discernible deformations for all three control techniques. Figures 12a 
and 12b show the total modal energy and actual command effort versus 
time, respectively, for three control implementation techniques. Figure 
12c illustrates the antenna rotation about the Xg direction for all 
three cases. The addition of just two actuators resulted in much better 
performance and a reduction in effort of approximately 50%. Hence, the 
maneuver and vibration control approach is successful even in extreme 
situations and the necessary number of actuators depends on the 
performance requirements. 


VIII. CONCLUSIONS 

The perturbation method of solving the equations of motion permits 
the maneuver and vibration controls to be formulated independently. 

This eliminates the numerical difficulties encountered Refs. 3-5 when 
many modes are controlled. The quasi-modal equations of motion permit 



an efficient vibration simulation as well as a straightforward vibration 
control formulation. Both uniform damping control and natural control 
were found to be effective in controlling vibration even during 
relatively quick maneuvers. The force distributions for both the 
maneuver control and uniform damping control depend exclusively on the 
mass matrix. Because the mass matrix is time invariant, these control 
formulations are not affected during the maneuver and, hence, are 
robust. The optimal actuator locations are the points of maximum mass 
both for maneuver and vibration controls. Actuator dynamics causes a 
time lag in force application, which degrades control performance. 
Including the actuator dynamics along with the structure dynamics in the 
control formulation minimizes this degradation. 
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Table 1 Natural Frequencies (Hz) 


Mode (Zero Gravity) 

0.95626111 

I. 02205468 
2.85798288 
4.12238565 
7.13573328 

II. 8067296 
14.4703039 
29.3765971 
31.8650183 

10 35.5681068 



Figure 1. Spacecraft in Orbit 
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Figure 3. Comparison of Maneuver Strategies 
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Figure 7. Comparison of Various Vibration Control Implementation 
Procedures for 30° Roll Maneuver 
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Figure 8. Antenna Rotation During 30° Maneuver 
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Figure 11. Effects of Actuator Saturation with Uniform Damping Control 
During 30° Maneuver 
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Figure 12. Implementation of 180° Maneuver with Various Numbers of 
Actuators 






